Minigap in superconductor— ferromagnet junctions with inhomogeneous magnetization 
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We consider the minigap in a disordered ferromagnet (F) in contact with a superconductor (S) in 
the situation when the magnetization of the F layer is inhomogeneous in space and noncolhnear. If 
the magnetization is strongly inhomogeneous, it effectively averages out, and the minigap survives 
up to the exchange field he ~ {L/a)ETh, where L is the thickness of the F layer, a is the scale on 
which the magnetization varies, and ETh is the Thouless energy. Technically, we use the "triplet" 
version of the Usadel equations, including both singlet and triplet components of the Green's func- 
tions. In many cases, the effect of disordered magnetization may be effectively included in the 
conventional Usadel equations as the spin-flip scattering term. In the case of low-dimensional mag- 
netic inhomogeneities (we consider spiral magnetization as an example), however, the full set of 
"triplet" equations must be solved. 

PACS numbers; 74.45.-|-c, 75.60.Ch, 74.78.Fk 



I. INTRODUCTION 

A normal metal in contact with a superconductor 
acquires some superconductive properties; this phe- 
nomenon is known as the proximity effect. One of the 
most prominent manifestations of the proximity effect is 
the (mini)gap in the single-electron spectrum of the nor- 
mal metal. While the density of states in the normal 
metal is spatially dependent, the minigap is a property 
of the normal part as a whole. In the diffusive limit, if the 
thickness of the normal metal L is larger than the coher- 
ence length, the characteristic scale of the minigap Eg is 
set by the Thouless energy E'xh ~ D/L^, where D is the 
diffusion constant (we put /i = Experimentally, the 
minigap can be directly probed by a scanning tunneling 
microscope (see, e.g., Ref. |2and references therein). 

Nowadays, the proximity effect in more complicated 
superconductor-ferromagnet (SF) structures is actively 
studied (see Ref. for a recent review of theoretical and 
experimental progress). While the majority of earlier 
theoretical studies treated the case of single-domain fer- 
romagnets, one of the open questions is the influence of a 
domain structure. For example, in some setups, inhomo- 
geneous magnetization may generate long-range triplet 
superconducting correlations 

In a single-domain ferromagnet, the exchange field h 
(measured in energy units) shifts the densities of states 
for the two spin subbands in the opposite directions, 
therefore the minigap in the spectrum closes at /i ~ Etyi^ 
A domain structure would effectively lead to averag- 
ing the nonuniform field h(r) acting on electrons; hence 
one can expect that the minigap would survive even at 
h ^ -Exh- In previous studies, the influence of inhomoge- 
neous magnetization on the density of states was consid- 
ered only in situations where the minigap was absent 
while the influence of a domain structure on the minigap 
has not been investigated. 

In this paper, we study the minigap in SF junctions in 
the presence of inhomogeneous magnetization with the 



help of the Usadel equations in the form allowing for the 
triplet superconducting component.^ We introduce a gen- 
eralization of the well-known 9 parametrization,® which 
involves, in addition to 9, a complex vector function M 
with the same number of components as for the exchange 
field h (three in the general case). 

The problem has three relevant energy scales: -Exh, h, 
and the energy Ed = D/a? associated with the length 
scale a over which the magnetization varies (typically, of 
the order of the domain size). Note that the supercon- 
ducting gap A will not play any role, as it is taken to 
be much larger than Exh and h. Throughout the pa- 
per we assume that the domains are small and that the 
ferromagnetic exchange field is weak. 



Ed > K 



Ed->E' 



Xh, 



(1) 



making no assumption on the relative scale of h and i?xh- 
First we consider a SF (or SFS) system with randomly 
oriented magnetization [Fig. ^a)] described by the pair 
correlation function 



h(r)h(r') = F{r - r'). 



(2) 



where the averaging is over an ensemble. We assume that 
the integral of the correlation function F(r — r') vanishes. 



J F{r) dr = 0, 



(3) 



and take a to be the corresponding length scale. Physi- 
cally, it means that the field h averages out on the scale 
a (of the order of several domain sizes). This assumption 
appears necessary for the smallness of the triplet super- 
conducting correlations and for the possibility to take 
them into account as an additional effective local term 
in the Usadel equation on the singlet component. Such 
a reduction is possible for the three-dimensional disorder 
(in the whole range of parameters, except at very low 
exchange fields and at energies near the minigap edge, 



2 



(a) 






s 


^\ I \ 


t 










L 





(c) 




u u u U U 



2L 




(d) 



FIG. 1: (a) SF junction with disordered magnetization, 
(b) Schematic dependence of the minigap Eg on the efTec- 
tive spin-fiip rate Fsf (see Ref. for details); the spin-flip 
approximation is not applicable in the shaded region (see the 
text for details), (c) and (d) SFS and SF junctions with spiral 
magnetic order. 



see discussion below): in this case, the inhomogeneous 
field h may be eff'ectively included in the conventional 
Usadel equation for the singlet component of the pairing 
correlations as the spin-flip term^ with the spin-flip rate 



(hV-2h) 



D 



Ed 



(4) 



(the averaging is over a region much larger than a) . Note 
that our result is different from that for pointlike impu- 
rities or disorderj2iiS because we consider h varying on 
length scales much larger than the elastic mean free path. 

The SF system with constant Tsf in the F part has 
already been studied in detail. The critical value of F^/ 
where the gap closes is Tsf/E^h = tt^/IO, which trans- 
lates into the critical value of the exchange field 



V ErhEd ^ — ETh- 



(5) 



If the disorder is one- or two-dimensional, one cannot 
easily reduce the nonhomogeneous field h to a spin-flip 
term and needs to solve the full nonlinear set of equa- 
tions. As an illustration we consider magnetic systems 
with the spiral magnetization [Figs.QIc) and^d)]. Re- 
markably, the critical values of the exchange fields for 
those systems also have the order of magnitude © . 



II. METHOD 

The Usadel equation for the Green's function g (which 
is a 4 X 4 matrix in the Nambu and spin spaces satisfy- 
ing the normalization condition = 1 = tq&q) has the 
forr^^^ 



{gVg) + iE [h^o,g] - i [r3(ho-), g] - A [fi^o. 



= 0, 
(6) 



where f and a are the Pauli matrices in the Nambu and 
spin spaces, respectively; E is the energy; h is the ex- 
change field; and A is the superconducting order param- 
eter. For simplicity, we choose A to be real and do not 
consider systems with a phase difference. 
The solution has the form 



g {go&o + go-) + fi (/oCTo + to-) 



(7) 



The normalization condition can be resolved by the 
parametrization 



go = Mq cos 0, 
fo = Mq sin 61, 



g = iMsin^, 

f = -iMcos6l, (8) 



with complex functions 6, Mq, and M, and with the con- 
straint 



1. 



(9) 



The Usadel equation © then yields one scalar and one 
vector equation: 



^^7^6* + A/o {iE sin 6* + A cos ( 



(hM) cos 61 = 0, 

(10) 



D 
"2" 



MV^Afo - A/oV^M) - M(i£:cos6' - A sin 6*) 
hA/osin6' = 0. (11) 



We may note the following general properties of Eqs. 
(Cni and (CH): 

(i) In the absence of ferromagnetism, h = 0, Mq = 1, 
M = 0, and Eqs. H10() and (|ll|l reduce to the conventional 
Usadel equation. 

(ii) For the uniform nonzero magnetization h = const, 
Eqs. ((Tn|l and ((TT|l imply the triplet vector M directed 
along the field h.^"^ 

(iii) A convenient feature of the parametrization (jSJ is 
that in the Matsubara representation [with —iE replaced 
by the Matsubara frequency ujn in Eqs. (|10|1 and l|ll|) ] the 
functions 9, Mq, and M are real. 

The density of states (summed over spin projections 
and normalized to the normal-metallic value) is stan- 
dardly expressed via the retarded and advanced Green's 
fimctions, f — Tr [f^ao [g^ — g^y\ /8, which yields 



Re go = Re(A/o cos 0) . 



(12) 



Below we consider SF and SFS systems and employ the 
rigid boundary conditions, which imply that the bulk so- 
lution in the superconductor with constant A is valid up 
to the SF interface, and the Green's function is continu- 
ous. The rigid boundary conditions are applicable if the 
SF interface is transparent and the F material is much 
more disordered than the S one. Thus at SF interfaces 



TT 

2' 



Afn 



1, 



M = 0, 



(13) 
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where the first condition is justified since /i <C A and we 
consider energies E < -Exhj hence E ^ A. The boundary 
conditions at the free surface of the F layer in the SF 
system are 



d9 

dz 



0, 



dMo 
dz 



0, 



dM 

dz 



0. 



(14) 



Below we consider Eqs. (IIUII and (|ll|l in the F part where 
A = 0; the superconducting correlations are induced due 
to the boundary conditions H13|l . 



III. DISORDERED MAGNETIZATION 

We start our analysis with the disordered system 
[Fig. ^a)]. To reduce the effect of the inhomogeneous 
field h to a spin-flip term, we use the self-consistent 
scheme for determining the spin-flip rate Tsj: We first 
solve the conventional "spin-flip" Usadel equation, 



D. 



-V^e + iE sin 6 - 2Vsf sin 6* cos 6* = 



(15) 



We expect that variations of h in space lead to the ef- 
fective averaging of the magnetization, hence the vector 
part M of the Green's function is small. Thus we substi- 
tute the solution of Eq. (|15|l into the linearized version 
of the vector Usadel equation (|ll|l to obtain 



D 



+iEcose M 



-hsin^. 



(16) 



This equation is solved for M with the appropriate 
boundary conditions, and further Tsf is found self- 
consistently as 



hM 

2 sin 6* 



(17) 



where the average is taken over a region much larger than 
the characteristic scale of variation of the field h. 

For this scheme to work, not only must we require 
|]V[| <C 1, but the (typically stronger) condition 



IMP < 



E' 



Th 



(18) 



[since in the scalar Usadel equation (|10|l we neglect the 
term quadratic in M, while keeping Tsf]. 

Equation H16|) is an inhomogeneous linear equation on 
M. We consider energies E < i?Th, hence \9\ ~ 1, and 
the {iE cos 9) term affects the Green's function of the 
linear operator in the left-hand side at a length scale of 
order L. It can be neglected if it produces a small relative 
correction to M, which in turn depends on the number 
of dimensions in which h(r) varies. While we always 
consider three-dimensional (3D) samples, h can either 
depend on all the three coordinates or be a function of 
only two (quasi-2D films) or one (quasi-lD wires) coordi- 
nate. Without the {iE cos 9) term, the Green's function 



of Eq. (Uni) is G ocl/r in 3D, G cx Inr in 2D, and G cx r 
in ID. The inaccuracy introduced due to neglecting the 
{iE cos 9) term is 



6M{r) = - j SG{r - ri)h(ri) sin6l(ri)d'^ri. 



(19) 



where the correction SG is negligible at length scales 
smaller than L. A straightforward estimate using the 
correlation function (0) with the vanishing integral Q 
gives ]JM|2 ~ a^+'^L^-'^K^/D^ and |Mp - a'^h^/D^, 
hence 



\SMl 
IMP 



(! 



(20) 



Therefore we can neglect the {iE cos 9) term in 3D, can 
obtain order-of-magnitudc estimates in this way in 2D, 
and cannot do it in ID. Note that we have also neglected 
the effect of boundaries. It can be taken into account via 
"mirror charges" (of the same or opposite sign for the 
boundaries at which dM/dx = or M = 0, respectively), 
and can be neglected under the same conditions. 

Now we assume that the disorder is three-dimensional 
and neglect the {iE cos 9) term in Eq. H16|l . Further, since 
h varies in space much faster than 9, we find 



M^--(V- 



^h) sin 61, 



(21) 



and Eq. H17() immediately leads to the effective spin-flip 
scattering rate (0J. That expression may be interpreted 
as the electrostatic energy of the "charges" h. Estimat- 
ing Tgf ^ h'^a? /D and |Mp ~ h^a^/D"^ translates our 
assumptions on the smallncss of M into the conditions 
(|T|l. which thus conflrm our linearization in M. 

However, to finally justify the spin-flip approximation, 
we must avoid one more danger: the noninvertibility of 
the "Hamiltonian" -[{D /2)\/'^ + iE cos 9]. The calcu- 
lations above were made without taking this possibility 
into account. At the same time, this happens exactly 
at the gap edge at Vgf = (i.e., at zero exchange field 
h). Indeed, below the gap 9 — 7r/2 -|- id, where i? is 
real. Parametrizing the solutions of the Usadel equation 
at Tsf = by ■i?05 the value of •& at the free F surface 
in the SF junction or in the center of the F layer in the 
SFS junction, and taking into account that the normal- 
metallic minigapi Eg = 0.78i?Th is the maximum value 
of i?(??o)j we see that dd/d^o is the zero-energy eigen- 
function of the above "Hamiltonian" . 

So in the vicinity of the point Tgf = 0, E — Eg, 
the spin-flip approximation breaks down. In this region, 
the zero mode becomes a low-energy mode with a small 
nonzero energy Eq. Expanding M = ^ qM^ in the ba- 
sis of the normalized eigenfunctions M^ of the "Hamilto- 
nian" on the left-hand side of Eq. (|16|l , we find that the 
contribution of the low-energy mode is determined by 



Co 



/ d'^r(hMo)sin6i. 

Eq J 



(22) 



Using the correlation function we estimate |cop ~ 
a^+'^L-'^h'^ /El. Finally, the condition |coMop < 
T sf / E'Yii that the contribution from the low-energy mode 
is small, yields 



'E^ ^ Vl/ 



(23) 



The energy Eq grows as we go away from the gap-closing 
point at Tsf = 0, E = Eg. The perturbation theory 
yields 



Eri 



\E-E„ 



E- 



Th 



E- 



Th 



E- 



Th 



(24) 



The combination of Eqs. (|23(l and l|24|l determines a small 
region near the T^f = 0, E — Eg point [shown as the 
shaded region in Fig.^b)], where the spin-flip approxi- 
mation (with the effective rate Tsf) breaks down due to 
the noninvertibility of the linear operator in Eq. (|16f) . 

Summarizing, the Usadel equations H1U|I and Ijlll) in 
the 3D problem reduce to the conventional Usadel equa- 
tion H15|) with the spin-flip rate Q if the conditions ^ 
and H23|) are satisfied. Then the results of Ref. 11 for 
the Eg{Tsf) dependence apply. The critical value of the 
exchange field, at which the minigap vanishes, is given 
by Eq. (|SJ) and is much larger than he in the case of a 
homogeneous ferromagnet. 



IV. SPIRAL MAGNETIZATION 

The above discussion in terms of the spin-flip approxi- 
mation applies to the 3D case and should give an order-of- 
magnitude estimate in the 2D case. In ID setups (with 
the field h depending only on one of the coordinates), 
one generally needs to solve the full nonlinear set of Eqs. 
(|10|l and Hll() . Below we present an analysis of the two 
examples [SFS and SF systems, see Figs, ^c) and^^d), 
respectively] with the regular spiral magnetic structure: 



h — h{cos kz, sin kz, 0), 



(25) 



where the z direction is perpendicular to the SF inter- 
face(s). In both the examples we are unable to ana- 
lytically find the full dependence of the minigap on the 
strength of the field h, but we calculate the critical value 
he at which the minigap closes. 

Technically, the SF case turns out to be more compli- 
cated, because the linearization over M does not work 
in this situation (while it is still applicable in the SFS 
system). 



A. SFS case 

We consider the spiral ferromagnet of length 2L 
[Fig.njc)] with the rigid boundary conditions (|13f) at the 
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FIG. 2: Effective quantum-mechanical potentials for (a) SFS 
and (b) SF junctions. The minigap vanishes when the Hamil- 
tonian with the corresponding potential has the ground state 
with zero energy. In both cases, the impenetrable walls are 
at 2; = and z = 2L. (a) Rectangular potential well, (b) Po- 
tential well with a 5-functional contribution in the middle. 



SF interfaces. With these boundary conditions and at 
zero energy it is possible to perform a derivation similar 
to that in the disordered 3D case. To find the gap-closing 
point, we find the minimal value of h at which the E = 
solution with 9 — 7r/2 has a bifurcation (a gapless so- 
lution forks out). Linearizing in (9 — Tr/2) and in M at 
E = 0, we solve Eq. 1)16(1 . without taking into account 
boundary conditions, by the oscillating function 



2h 



(26) 



The boundary conditions M(0) = M(2L) = are easily 
satisfied by adding a smooth (nonoscillating) term linear 
in z of the same order of magnitude: M = Mqsc + ^sm- 
On substituting this solution into the linearized scalar 
Usadel equation, we find the effective equation for the 
smooth part of 9{z) [the osciUating part of 9{z) has a 
smaller order of magnitude and is neglected; only Mqsc 
needs to be taken into account at this step]: 



= 



with the boundary conditions 

61(0) = 9{2L) = 



(27) 



(28) 



Equation H27|l can be viewed as the Schrodinger equation 
at zero energy with h ~ I, mass 1/2, and the constant 
potential -Ah"^ /{D'^k^), while {9 - 7r/2) plays the role of 
the eigenfunction. The boundary conditions H28() corre- 
spond to the impenetrable walls at z — and z — 2L, 
see Fig.|2i;a). 

Therefore, the bifurcation of a nontrivial solution for 
the problem H27() and 1(281) (i.e., closing of the minigap) 
with increasing h corresponds to the ground-state energy 
crossing zero. This immediately leads to the result 



lSFS 



7r\ Dk 



L 



(29) 
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We verify that the hnearization over M is justified under 
the conditions (Q) (with ~ Dk^) which, in turn, re- 
duce to the single requirement fcL 3> 1 at /i . Note 
that Eq. (|29|l formally coincides with the expression (HI) 
in combination with the critical value of Tsf from Ref.illi 
The critical value of the exchange field is, as in the pre- 
vious example, much larger than the Thouless energy: 
/if s ^ (fcL)iJTh. 



B. SF case 

Another example we consider is one half of the above 
SFS junction: the SF junction with the ferromagnet of 
length L [Fig.^d)]. In this case, the open boundary con- 
ditions cannot be satisfied at small M and the lineariza- 
tion over M does not work. [Physically, M is enhanced 
because electrons incident on the outer F surface through 
the field h(z), feel the same h(z) after reflection, hence 
the effective average of h is nonzero.] Thus we need to 
solve the full nonlinear Eq. Ullfl aX E = 0. Below the 
gap, M is real, and it is convenient to introduce the new 
complex function 



m{z) = Ml -I- iM2 



(30) 



(where Mi and M2 are the two components of the vector 
M). Then in the new complex notation Eq. takes 
the form 



he''''' = 0, Mo = VT 



D 
"2 

(31) 

with the boundary conditions m(0) = and m'(L) = 0. 
We separate rapidly oscillating modes, 



^bo + hie' 



-le 



— ikz 



(32) 



and keep only the leading ones (the amplitudes 6„ are 
slow functions of z and decrease in magnitude with in- 
creasing |n|, as we verify below). The equations for differ- 
ent modes are coupled, and after a straightforward alge- 
bra, we express the amplitudes 61 and &_i via the smooth 
part hQ-. 



hi 



Dk^ 



(2 



\h,?) 



h-i = 



—h^ 



(33) 



Substituting these amplitudes in Eq. I|31|l and 



parametrizing 60 = 



ikL 



sinh ip, we arrive at the 



■ sinh </? cosh = 0, 



(34) 



effective equation on f{z): 

with the boundary conditions 

(p(0) = 0, ip'{L) = X cosh ip{L), (35) 

where k = 2h/{Dk). The solution to this equation must 
be substituted back into the scalar Usadel equation H1(J|) 



to find the bifurcation point of the solution 9 = tt/2. 
Linearizing the equation in the vicinity of the bifurcation 
point, we obtain 



sinh^ if 



(/or sinh (p) sin k{z — L) 



o sinh ifi^, 
^ 1 cos 2k{z - L) + 



(36) 



Again, we need to take into account not only the smooth 
part of 6, but also the oscillating modes at ±fc: 



TT 

2 



to + ti sin k{z — L) + 



(37) 



Solving coupled equations for the amplitudes i„ (which 
are slow functions of z and decrease in magnitude with 
increasing the order of the harmonic), we find, to the 
leading order in >?/fc. 



(|sinh<,) 



to 



(38) 



and finally obtain the equation for the smooth component 
to- 

{Vl + cosh^ ip) to = 0, (39) 
with the boundary conditions 

to(0) = 0, t'aiV) - [^sinh(^(L)]io(i). (40) 

At the bifurcation point (i.e., when the minigap closes), 
this homogeneous linear equation acquires a nonzero so- 
lution. The resulting problem again has a quantum- 
mechanical analogy if we symmetrically continue the po- 
tential cosh^ 9?(z) in the Schrodinger equation l|39|) 
from the (0, L) interval to (L,2L); see Fig. |2Ib). As a 
result, we obtain the Schrodinger equation at zero en- 
ergy with h—X and mass 1/2, while to plays the role of 
the eigenfunction. At z = and z — 2L, we obtain the 
impenetrable walls, while the boundary condition (|4()|l at 
z = L corresponds to the following (5-functional contri- 
bution to the potential: —2 [>irsinhiy9(L)] 5{x — L). 

Therefore, the appearance of a nontrivial solution for 
the problem (|39|l and H40|) at the bifurcation point with 
increasing h corresponds to the situation, in which the 
ground state of the Schrodinger equation has zero energy. 
Numerically, we find that this happens at = 0.5955/L, 
which translates into the critical value of the exchange 
field 



/if = 0.2977:?^. 



(41) 



Having found /if, we can check that 60 ^ 1, &±i ^ 1/kL, 
h±2 ~ l/(/eL)^, ti ^ to/kL, etc., confirming the validity 
of the expansions H32() and (|37|l . Note that the critical 
strength of h has the same order of magnitude as in the 
SFS case [and thus also agrees with the estimate ©J, but 
it has a smaller numerical prefactor. 
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V. CONCLUSIONS 

In conclusion, we have shown that SF and SFS systems 
with random three-dimensional domain disorder admit 
an effective description in terms of the spin- flip rate Tgf ■ 
We have also considered SF and SFS junctions with spiral 
magnetic order, where the spin-flip approximation does 
not hold. In all the systems considered, the minigap sur- 
vives up to the exchange fields of the order (|3J), i.e., much 
larger than E'Th • For an experimental observation of the 
effects described in the present paper, it is crucial to use 
weak ferromagnets with a very fine domain structure, so 



that the main constraints 1^ are satisfied. 
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